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Abstract 

The transverse momentum distribution of single vector boson produc- 
tion at hadron colliders provides useful ways of testing the Standard Model 
and searching new physics beyond the Standard Model. We study large p T 
hadroproduction of Z-boson as a probe of gluon distributions inside proton. 
We investigate how to get initial gluon-involving contributions, or how to 
subtract quark-quark (or -antiquark) contributions from total cross section. 
We also investigated the simultaneous measurement of the rapidity and the 
transverse momentum of the produced Z boson, to obtain momentum frac- 
tions of initial partons. And we extracted relevant uncertainties involving in 
experimental and theoretical analyses. This large p T hadroproduction of Z 
can be used as constraints on analyses of global parton (gluon and quarks) 
distribution functions inside proton. 
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1 Introduction 



The transverse momentum distribution of single vector boson production provides 
useful methods of testing the Standard Model and searching new physics beyond the 
Standard Model. Here we would like to investigate mainly the gluon distribution 
inside proton using large p T production of Z-boson at hadron colliders. 

The recent improvements in the data of deep-inelastic lepton-nucleus scattering 
and of the Drell-Yan process have allowed a definitive determination of the quark 
distribution functions inside the nucleon. However, the gluon distribution G(x, Q 2 ) 
is not well constrained by these processes, since it only enters as a second-order 
effect. On the other hand, the gluon contributes to the lowest order for vector 
boson hadroproduction with large transverse momentum and it has become normal 
practice to use fixed-target pp(p) — > VX data to determine the gluon distribution for 
x ~ 2p T / y/s (where V represents a vector boson and p T is the transverse momentum 
of the produced vector boson). The prompt photon production pp — > 7X has 
been already used on that purpose 0. In principle, collider data obtained for 
pp — ¥ 7X at the Fermilab energy (y/s = 1.8 TeV) with p T ~ 10 GeV could probe 
the gluon at x ~ 10~ 2 . However, a detailed study of uncertainties in the theorectical 
predictions of prompt-photon production at collider energies has shown that such 
determination will be difficult for some reasons p, |J . One problem is the importance 
of the bremsstrahlung component at small p T , in which the photon is radiated 
from an outgoing quark and so occurs in the debris of a hadronic jet. A second 
ambiguity is related to the choice of scales in the parton distribution functions and 
the QCD coupling strength a s . For these uncertainties, massive vector boson W 
or Z hadroproduction is more useful than 7. The bremsstrahlung W, Z processes 
are almost negligible; that is the outgoing quark (or anti-quark) will very rarely 
fragment into W or Z. And the choice of scales can be fixed more acceptably, for 
instance, as Q 2 = or M§. 
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Previously, as an example, it was pointed out that W production could be used 
as measures for gluon and heavy quark distributions |4[], which concentrated on the 
following subprocesses; 

(i) sg — > cW 
iii) eg — > cZ. 

They proposed (i) as a measure of the gluon, and the ratio of (ii) to (i) to determine 
the ratio c/s of quark densities. However, a W^-boson possibly decays leptonically 
into eu. And because one cannot detect neutrinos, one could confirm W production 
only through detecting an electron and large missing p T , which usually results in 
large systematic uncertainties. On the other hand, a Z-boson decays leptonically 
to e + e~ , so we can confirm Z production with less systematic uncertainties. Pacing 
with considerable new data of Z hadropro duct ion at Tevatron, we would like to 
investigate gluon densities of proton by using large p T hadroproduction of Z. At 
a glance, it seems possible to use the above subprocess, (ii) eg — > cZ, as a probe 
of gluon. But it is not so effective, since this subprocess involves the initial charm 
quark, and the resulting cross section would be very small. As one can see more 
details later, we instead use initial qg subprocesses, where q represents any of three 
light quarks. 

In this letter, we present the transverse momentum and rapidity distributions of 
Z boson, calculated up to the second order in QCD coupling. When the produced 
Z has small transverse momentum p T , there are large logarithms log(Q 2 j 'p 2 T ) , and 
perturbation theory for da/dp 2 T breaks down there. And one must either perform 
a resummation or restrict attention to the total cross section by integrating ana- 
lytically in p T . These techniques have been carried out to the first order in QCD 
||, and part of the second-order terms have been analyzed in the p T — > limit 0. 
Therefore, we here restrict ourselves to the region of large p T to avoid these difficul- 
ties in small p T . In the references |7], ||, the full analytic formulae of da / dp 2 T dy for 
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large p T vector boson production up to the second order in QCD had been already 
presented. Based on those formulae we study the gluon distribution inside proton 
through large p T hadroproduction of Z-boson. By including the full aa 2 contribu- 
tions, we are able to considerably reduce the theoretical errors, normally associated 
with the leading order (0(aa s )) results. We calculated K = da-iiZ) j da\{Z) as a 
function of p T (Z\ where da n denotes the differential cross section da/dp T (Z) in- 
cluding all QCD subprocesses up to 0(aa") at \/s =1.8 TeV using Q 2 = Mf , and 
found that K ~ 1.3 ~ 1.4 for p T (Z) > lOGeV. The corrections are significant, but 
show that the perturbative expansion is reasonable. 

Up to the order of aa 2 , pp — > ZX process includes initial qq (qq or qq), qg (qg) 
and gg scatterings in partonic level. Among these subprocesses, only qg (qg) and 
gg scatterings are initial gluon-involving parts. Therefore, one can analyze the 
gluon distribution by subtracting initial quark-quark contributions from the total 
subprocesses of Z hadroproduction. At the Tevatron energies the contribution from 
initial gg part is very small, and we can ignore in the numerical analyses. 



2 Detailed investigations of gluon distribution 

As mentioned in Section 1, we would like to use the formulae of reference for our 
detailed analyses. The authors of J7| presented the hadronic cross sections up to the 
order of aa 2 for A + B — » 7* + jet(s), 

do ' AB x — ^ f , , „j, „o, .r, „o> sdu 



J dx 1 dx 2 f t A (x 1 ,Q 2 )ff(x 2 ,Q 2 )^(a s (Q 2 ),x 1 P 1 ,x 2 P 2 ) 



dp 2 T dy 

where A and B are the initial hadrons, Pi are their momenta, and da is the factorized 



partonic cross section in the MS factorization scheme. 

The large p T hadroproduction of Z up to aa 2 can be calculated by the charac- 
teristic five groups of diagrams, as follows 

(i) qq -> gZ, 



(ii) Next-to-leading diagrams for qq — > gZ, 
(Hi) qq -> ggZ, 

(iv) qq — > qqZ, 

(v) qq — > qqZ. 

Diagrams (i) are the first-order ones for Z production through qq annihilation. Parts 
of the second-order contribution come from the interference of these diagrams with 
the one-loop corrections of diagrams (ii). The rest comes from 2- jet productions 
accompanying the Z, such as the diagrams (Hi). The diagrams through initial qg- 
scattering can be obtained from (i), (ii) and (Hi) by crossing them. The diagrams 
through gg-scattering can be similarly obtained from (Hi) alone. The diagrams 
(iv) give the remaining parts to gg-scattering. And the diagrams (v) are for qq- 
scattering for Z hadroproduction. Therefore, we have three types of contributions 
for Z hadroproduction 

• QQ — > Z + 1 or 2 jets, 

• QG -> Z + 1 or 2jets, 

• GG -> Z + 2jets, 

where Q represents initial parton of a quark or anti-quark, and G for a gluon. 

To further investigate our observations, we numerically calculate the transverse 
momentum p T , and the rapidity y distributions of Z at the Tevatron energies in- 
cluding all QCD subprocesses up to the order aa 2 s in pp collisions. We use three 
input parton distribution functions; MRS (A) §, GRV(HO) |IJ and CTEQ3 Jll|. 
Fig. 1 shows p T distributions using three parton distributions at yfs = 1.8 TeV, 
where we use Q 2 = Mf , and include three light quarks (u, d, s) as initial quark- 
partons. As mentioned earlier, initial GG-scattering part is very small, and we do 
not show it explicitly here. (Of cause, we include GG-part to calculate total cross 
section.) We find that their gluon distributions are almost same at Q 2 ~ Mf , and 
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in Fig. 1 the deviations of distributions are mainly due to the difference of sea-quark 
distributions. We also note from Fig. 1 that it looks very important to subtract 
the initial QQ-contributions to extract contributions involving gluons only. How- 
ever, the pp — > Z + jets process contains many different subprocesses, and it might 
be impossible to distinguish the initial QQ-contributions from the QG-oiaes exper- 
imentally by analyzing shapes of final hadronic jets, and etc. Now that we have 
theoretical predictions of the initial QQ-contributions, we can in principle subtract 
those from the experimental total result. And by comparing the remaining experi- 
mental cross section and theoretical one, one can probe the gluon distributions inside 
proton. However, as can be seen in Fig. 1 for p T > 10 GeV, QQ-contributions are 
about 70% of total cross section, and it looks inappropriate to subtract QQ-p&rt, 
which is not a minor part, but a major contribution. Therefore, one should raise 
p T minimum cut up to about 40 ~ 50 GeV, to make QQ- and QG-contributions at 
least experimentally comparable. These practical considerations will be presented 
in more detail in Section 3. 

Assuming that we could experimentally distinguish Z + ljet from Z + 2jets, let 
us first consider simple Z + ljet case. In partonic level, pp — > Z + ljet process 
includes only two types of diagrams; 

• qq -> Zg, 

• qg -> Zq (or qg -> Zq). 

In these processes, we can extract initial gluon-involving contributions by subtract- 
ing theoretical prediction of qq — > Zg subprocess from the experimental total results 
of pp — > Z + ljet. Futhermore, if we could distinguish the quark-jet from the gluon- 
jet through theoretical and experimental combined analyses, only the initial gluon 
contributing part can be extracted experimentally from Z + 1 quark-jet. Fig. 2 
shows the p T distributions of Z + ljet from qq — > Zg and qg — ■> Zq (qg — > Zq). For 
comparison, we also show p T distributions of pp — > Z + ljet and Z + 2jets. 
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Until now we have only considered the transverse momentum p T distributions. 
We now consider analyses with simultaneous measurement of both rapidity y and 
p T of the produced Z boson. The y and p T of the Z boson in the laboratory frame 
are related to the momentum fractions of initial partons by 
1, E + p L 

y = o m 



2 E-p 



L 

2 



, Mf + xix 2 s - s 2 (Mf + xix 2 s - s 2 ) 2 

where s is the invariant mass of the incoming hadrons, and s 2 is the invariant mass 
of final two-jets when two jets accompany the produced Z. And xi, x 2 are the four- 
momentum fractions of the colliding partons. If we consider only Z + ljet processes, 
the value of s 2 in the above relation becomes zero. By using the above relation we 
can directly obtain the gluon momentum fraction inside incoming proton. 



3 Discussions 

We now study relevant uncertainties in analyzing large p T hadropro duct ion of Z, 
which one must take into account in theoretical and experimental analyses. The 
first uncertainty comes from QCD-scale dependences. Although there is much less 
ambiguity associated with the choice of scales compared to direct photon 
explained before, there still remains weak QCD-scale dependence because of con- 
sidering only finite order of perturbative corrections. In Table 1 we give total cross 
sections for pp — > ZX production integrated over the region p T (Z) > p^ in = 10 
GeV for various Q 2 . The value of p^in ls chosen so as to retain as many clear events 
as possible and yet to remain in a region where perturbation theory gives a reliable 
prediction. The same scale is chosen for the parton densities and QCD strength 
ct s (Q 2 )- But these scales need not be the same. Table 2 shows the sensitivity, when 
we choose different scales for a s (Q 2 ) and parton distribution functions. We note 
that there show significant differences in results between Q 2 = Mf and Q 2 = p 2 T /2 
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Table 1: The integrated total cross sections with the same scales for a s and parton 
distributions. 



a (nb) 



p T >lOGeV, \y\<2.5 



scale Q 2 


M 2 


p\n 


P 2 T 


2 P 2 T 


MRS (A) 


1.84 


2.16 


2.06 


1.99 


GRV(HO) 


1.73 


2.05 


1.95 


1.88 


CTEQ3 


1.87 


2.19 


2.10 


2.02 



for parton distributions, where we fix the scale of a s (Q 2 ) as Mf, as in Table 2. 
Therefore, once we have sufficient experimental data for this process (pp — > ZX), 
we could even investigate scale dependence of QCD-coupling a s , as well as parton 
distribution functions. 

Table 2: The integrated total cross sections with different scales for a s and parton 
distributions. We fix Q 2 = Mf for a s and use four different scales for parton 
distribution functions. 



a (nb) 



p T >10GeV, \y\<2.5 



scale Q 2 


Mf 


pin 


pI 


2p 2 T 


MRS (A) 


1.84 


1.52 


1.59 


1.64 


GRV(HO) 


1.73 


1.43 


1.50 


1.54 


CTEQ3 


1.87 


1.54 


1.61 


1.66 



Next we still have uncertainties related to structure functions of initial quarks, 
when we try to derive informations on gluon. Here we use three different parametriza- 
tions; MRS(A), GRV(HO), CTEQ3. Martin et. al. [|§ improved their parametriza- 
tion of parton densities through a new global analysis on deep inelastic scatterings 
and related data including the recent measurements of F2 at HERA, on the asym- 
metry of the rapidity distributions of production [|1^] at Tevatron, and on the 
asymmetry in Drell-Yan production in pp and pn collisions ||13|1 . We use their latest 



version of MRS (A) functions. CTEQ Collaboration [11] also improved their distri- 
bution functions, incorporating several new types of data and presented the new 



version CTEQ3, which we use here. Gliick et. al. [n| predicted the parton distri- 
butions down to x ~ 10~ 4 and Q 2 ~ 0.3 GeV 2 , using the data from deep inelastic 
scattering experiments at x > 10~ 2 together with the idea that at some low reso- 
lution scale the nucleon consists entirely of valence quarks and valence-like gluons. 
They presented parton distributions obtained for the leading order (GRVLO) as well 
as for the higher order (GRVHO) calculation. We use GRVHO. 

The simultaneous measurement of y and p T would enable us to determine the 
momentum fractions of initial partons from Eq. (1), and it was previously pointed 
out that one should be very careful in dealing with s 2 , where s 2 is the invariant mass 
of final state two jets, 

$2 = (Pjetl+Pjet2) 2 - 

If only one jet presents in the final states, s 2 = in our massless parton level approx- 
imation, as explained before. Experimenalists should find here the best optimized 
solution through realistic Monte Carlo studies to differentiate one-jet from two-jets 
final states, which are accompanying the produced Z-boson. Fig. 3 shows the p T 
distributions of pp —>■ Z + jets for y = 0. 

We next calculate numerically the number of events at Tevatron energies. For 
\y\ < 2.5 and p T > 10 GeV, we get 

Number of events /lOpb' 1 = [2.00 ± 0.15] x 10 4 for MRS (A), 

= [1.89 ± 0.16] x 10 4 for GRV(HO), 
= [2.03 ± 0.16] x 10 4 for CTEQ3, 

where the errors are such that the lower and upper limits correspond to the scale 
choice of Q 2 = M§ and Q 2 = p 2 /2, respectively. As previously mentioned, about 
70% of the above events is from the initial QQ-contributions, and QG-contributions, 
which we are interested in for our investigation of gluon distributions, are only 30% 
of total events. It would not be so effective to subtract from total number of events 
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QQ-contributions, which are much larger than QG-ones. Therefore, we need raise 
p T minimum cut to make initial QG-contributions at least comparable to QQ-ones. 
Fig. 4 shows the rapidity distributions with (a) p T > 10 GeV, and (b) p T > 40 GeV. 
And for \y\ < 2.5 and p T > 40 GeV, we get 

Number of events /10pb- 1 = [2.27 ± 0.07] x 10 3 for MRS (A), 

= [2.18 ± 0.07] x 10 3 for GRV(HO), 
= [2.35 ± 0.08] x 10 3 for CTEQ3. 

As we can see in Fig. 4(b), about a half of the total number of events is from 
initial QG-contributions. Therefore, we can summarize, assuming luminosity C = 
lOpb' 1 /year- 
Number of events/year (p T > lOGeV) = [1.96 ± 0.08 ± 0.16] x 10 4 , 
Number of events/year (p T > AOGeV) = [2.27 ± 0.08 ± 0.08] x 10 3 , 

where the first error is due to structure functions, and the second is from scale (Q 2 ) 
dependence. For p T > 10 GeV (see Fig. 4(a)), the total errors are about 9%, and the 
QG-contributions are only about 30%. One the other hand, we note that for p T > 40 
GeV the initial QG-contributions are comparable to the QQ-ones (see Fig. 4(b)), 
and the theoretical errors are also reduced to 5%. Therefore, for large p T minimum 
cut we can investigate gluon structure functions with less theoretical uncertainties. 
Moreover, we can see from Fig. 4(b) that QG-contributions are even larger than QQ- 
ones at large rapidity region. It means that (Q(x q ) + G(x 9 ))-scattering contributions 
are dominant at large (x q — x g ), where x q and x g are momentum fractions of intital 
scattering quark and gluon, with the kinematic constraint of sx q x g > Mf . 

Fig. 5 shows rapidity distributions at y/s = 14 TeV (LHC energies). At these 
energies, QG-contributions exceed 80% of total cross section. If we consider only 
Z + ljet case at these LHC energies, because initial QG-scatterings produce only 
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quark-jet final states, we can investigate quark fragmentation functions using this 
large p T hadroproduction of Z-boson. For p T > 50 GeV, total number of events at 
= 14 TeV is 

Number of events /lOOpb' 1 = [3.86 ± 0.06 ± 0.11] x 10 5 , 

where the first (second) errors are due to structure function (QCD-scale) depen- 
dence. 

This process has been already used for other purposes. In recent papers jn| , it is 
reported that the cross sections for W and Z production in pp collisions at a/s =1.8 
TeV are measured at the Fermilab Tevatron colliders for final states; W — > ev e , 
Z — > e + e~, W — > [xv and Z — > And it is pointed out that assuming the 

Standard Model couplings, this result can be used to determine the width and mass 
of the W boson. 
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Figure Captions 



Fig.l The p T distributions of pp -> Z + X at ^/s = 1.8 TeV using Q 2 = Mf . QQ 
and QG represent the initial quark — quark (or —antiquark) and quark— gluon 
contributions respectively, and total represents for the sum of all contributions. 
The solid line is for MRS (A), the dash-dotted is for CTEQ3, and the dashed 
is for GRV(HO). 

Fig. 2 The p T distributions of pp — > Z + I jet from (QQ) qq — > Zg and (QG) 
qg — > Zc/ (55 — > Z<?) at y^s = 1.8 TeV using Q 2 = Mf . Also shown are the p T 
distributions of pp — * Z + ljet and Z + 2jets. 

Fig. 3 The p T distributions of pp — > Z + X for y = 0, at — 1-8 TeV using 
Q 2 = Mf. Also shown are the p T distributions of QG — > Z + ljet and 
Z + 2jets cases for y = 0. 

Fig.4 The rapidity distributions of QQ -> Z + X and QG -> Z + X (a) forp T > 10 
GeV, and (b) for p T > 40 GeV, at y/s = 1.8 TeV using Q 2 = Mf . 

Fig. 5 The rapidity distributions of — > Z + X and QG — > Z + X for p T > 50 
GeV at v/i = 14 TeV using Q 2 = Mf . 
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